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Abstract 

Condensed matter systems, such as acoustics in flowing fluids, light 
in moving dielectrics, or quasiparticles in a moving superfluid, can be 
used to mimic aspects of general relativity. More precisely these sys- 
tems (and others) provide experimentally accessible models of curved- 
space quantum field theory. As such they mimic kinematic aspects of 
general relativity, though typically they do not mimic the dynamics. 
Although these analogue models are thereby limited in their ability 
to duplicate all the effects of Einstein gravity they nevertheless are 
extremely important — they provide black hole analogues (some of 
which have already been seen experimentally) and lead to tests of ba- 
sic principles of curved-space quantum field theory. Currently these 
tests are still in the realm of gedanken- experiments, but there are plau- 
sible candidate models that should lead to laboratory experiments in 
the not too distant future. 
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1 Introduction 

Analogue models of (and to some extent for) general relativity have recently 
become a growth industry |[T|. Typically based on various condensed- matter 
systems, these analogue models are most often used for devising gedanken- 
experiments that probe the structure of curved-space quantum field theory. 
More boldly, they seem promising routes to providing real laboratory tests of 
the foundations of curved-space quantum field theory. (The most spectacular 
suggestion along these lines is that analogue models may make experimental 
tests of the Hawking radiation phenomenon a realistic possibility.) 

Ideas along these lines have, to some extent, been quietly in circulation 
almost since the inception of general relativity itself. Walther Gordon (of the 
Klein-Gordon equation) introduced a notion of "effective metric" to describe 
the effect of a refractive index on the propagation of light 0. The Russian 
school, as epitomized by Landau and Lifshitz, used notions developed in 
optics to represent gravitational fields in terms of an "equivalent refractive 
index" There is an extensive, but largely neglected samizdat literature 
(of extremely variable quality) that explores these issues. (For an extensive, 
though still not comprehensive, bibliography see @].) 

The modern revival is due largely to Unruh (and to some extent Mon- 
crief 0) who in the early eighties considered the use of hydrodynamic ana- 
logues, in which sound waves in a flowing fluid are mapped into a suitable 
scalar field theory in an effective curved spacetime — the "acoustic geome- 
try". (The precise statement, as will be described more fully below, is that 
sound in an irrotational inviscid barotropic fluid is identical to a massless 
minimally coupled scalar field in curved spacetime; and quantized sound [the 
phonon field] is identical to curved-space quantum field theory.) 

The nineties saw considerable work on the nature of Hawking radiation in 
these analogue models, still largely with the attitude that one was performing 
gedanken- experiments. It is only now, at the turn of the millennium, that 
serious consideration is being given to the actual construction of laboratory 
experiments. Three classes of system stand out as being the most likely to 
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lead to useful experimental probes: 

• Acoustics in Bose-Einstein condensates. 

• Slow light. 

• Quasiparticles in superfluids. 



2 Acoustics in BECs 

In this mini-survey we will mainly concentrate on acoustics in BECs, and give 
some feel for where we stand and what the near-term prospects are. Acoustic 
analogues of black holes are formed by supersonic fluid flow |^ . The flow 
entrains sound waves and forms a trapped region from which sound cannot 
escape. The surface of no return, the acoustic horizon, is qualitatively very 
similar to the event horizon of a general relativity black hole; in particular 
Hawking radiation (in this thermal bath of phonons with temperature 

proportional to the "surface gravity") is expected to occur 0. There are 
at least three physical situations in which acoustic horizons are known to 
occur: Bondi-Hoyle accretion the Parker wind (coronal outflow from 
a star), and supersonic wind tunnels. Recent improvements in the creation 



and control of Bose-Einstein condensates (see e.g., pO| , |TT[) have lead to 
a growing interest in these systems as experimental realizations of acoustic 
analogs of event horizons. In reference |]12| we considered supersonic flow 



of a BEC through a Laval nozzle (converging-diverging nozzle) in a quasi- 
one-dimensional approximation. We showed that this geometry allows the 
existence of a fluid flow with acoustic horizons without requiring any special 
external potential, and we then studied this flow with a view to finding situ- 
ations in which the Hawking effect is large. We were able to present simple 
physical estimates for the "surface gravity" and Hawking temperature, and 
so to identify an experimentally plausible configuration with a Hawking tem- 
perature of order 70 n K; this figure should be contrasted with the critical 
condensation temperature which is of the order of 90 n K. We stress that 
in present day experiments the actual physical temperature of the conden- 
sate, although difficult to measure, is believed to he well below this critical 
temperature. 
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2.1 From Gross— Pitaevskii to hydrodynamics 

Bose-Einstein condensates are most usefully described by the nonlinear Schro- 
dinger equation, also called the Gross-Pitaevskii equation, or sometimes the 
time-dependent Landau-Ginsburg equation: 

- iTi d^t, x) = -—V^^/j{t, x) + X\ iP{t, x). (2.1) 
im 

(We have suppressed the externally applied trapping potential for algebraic 
simplicity. For many technical details, and various extensions of the model, 
see Barcelo et al. [|T3l. That reference also contains an extensive back- 



ground bibliography.) Now use the Madelung representation to put the 
Schrodinger equation in "hydrodynamic" form: 

ip = ^/p exp{-i9 m/h). (2.2) 

Take real and imaginary parts: The imaginary part is a continuity equation 
for an irrotational fluid flow of velocity v = V6 and density p; while the real 
part is a Hamilton-Jacobi equation (Bernoulli equation; its gradient leads to 
the Euler equation). Specifically: 

dtp + V- {p Ve) = 0. (2.3) 

That is, the nonlinear Schrodinger equation is completely equivalent to irro- 
tational inviscid hydrodynamics with a particular form for the enthalpy 

h=['-l^^, (2.5) 

J p m 

plus a peculiar derivative self-interaction: 

The equation of state for this "quantum fluid" is calculated from the enthalpy 

The corresponding speed of sound is 

= = ^. (2.8) 

dp m 
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2.2 Acoustic metric 



To now extract a Lorentzian geometry, linearize around some background. 
In the low-momentum limit it is safe to neglect Vq. It is a by now standard 
result that the phonon is a massless minimally-coupled scalar that satisfies 
the d'Alembertian equation in the effective (inverse) metric 0, [l^ 



Po 



Here 



2 Po 

c„ = : 

" m ' 



(C^ I-VO0 Vq) 



(2.9) 



(2.10) 



It cannot be overemphasized that low-momentum phonon physics is com- 
pletely equivalent to (scalar) quantum field theory in curved spacetime. That 
is, everything that theorists have learned about curved space QFT can be 
carried over to this acoustic system, and conversely acoustic experiments 
can in principle be used to experimentally investigate curved space QFT. In 
particular, it is expected that acoustic black holes (dumb holes) will form 
when the condensate fiow goes supersonic, and that they will emit a thermal 
bath of Hawking radiation at a temperature related to the physical acceler- 
ation of the condensate as it crosses the acoustic horizon |^, |1^, |15[. For 
completeness we mention that the metric is 



X) 



Po 



- c: 



-Vq 



-Vo 
I 



(2.11) 



so the space-time interval can be written Q 



~cl dt^ + \ \dx-vo dt\ 



(2.12) 



The low-momentum phonon physics looks completely Lorentz invariant. (This 
is an acoustic Lorentz invariance mind you, with the speed of sound doing 
duty for the speed of light M.) 



^ Although Eq. 2.12 seems to imply that a standard BEC can only simulate metrics 
with conformally flat spatial sections, it can nevertheless be shown that if the condensate 
is characterized by some anisotropic mass tensor (realized, e.g., via some doping gradient) 
then non-conformally flat spatial sections could also be simulated . 
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2.3 Hawking effect 



2.3.1 Laval nozzle 

A general problem with the experimental construction of acoustic horizons 
is that many of the background fluid flows so far studied seem to require very 
special fine-tuned forms for the external potential. (See e.g., the Schwarzschild- 
like geometry in reference ||16|.) In this respect a possible improvement to- 
ward the realizability of acoustic horizons is the construction of a fiow in a 
trap which "geometrically constrains" the flow in such a way as to replace 
the need for a special external potential. An example of such a geometry 
is the so called Laval nozzle (converging- diverging nozzle). In particular we 
shall consider a pair of Laval nozzles; this provides a system which includes 
a region of supersonic flow bounded between two subsonic regions. Consider 



Figure 1: A pair of Laval nozzles: The second constriction is used to bring 
the Huid How back to subsonic velocities. 

such a nozzle pointing along the z axis. Let the cross sectional area be de- 
noted A{z). We apply, with appropriate modiflcations and simpliflcations, 



the calculations of references |T2| and |T6[. The crucial approximation is 
that transverse velocities (in the x and y directions) are small with respect 
to velocity along the z axis. Then, assuming steady flow, we can write the 
continuity equation in the form 

p{z) A{z) v{z) = J; J = constant. (2-13) 

The Euler equation (which we simplify by excluding external forces d$/dz, 
and excluding internal viscous friction /„) reduces to 

pv — = --f. (2.14) 

dz dz 
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Finally, we assume a barotropic equation of state p — p{p) , and define X' — 
dX/dz. Then continuity implies 



P ^-P 



while Euler implies 



{AvY 
(Av) 



-P 



A' v'' 
A V 



-P 



A' a 
A 



dp , 
pa = P- 
dp 



(2.15) 



(2.16) 



Defining the speed of sound by (? = dp/ dp, and eliminating p' between these 
two equations yields a form of the well-known "nozzle equation" 



2 2 

(,2 _ y2 



A' 



(2.17) 



The presence of the factor — v'^ in the denominator is crucial and leads to 
several interesting physical effects. For instance, if the physical acceleration 
is to be finite at the acoustic horizon, we need 



A' 



0. 



(2.18) 



This is a fine-tuning condition that forces the acoustic horizon (technically, 
the acoustic ergosurface) to form at exactly the narrowest part of the nozzle. 
(If external body forces and internal friction are not neglected, then there is 
a precise relationship between these forces and the location of the horizon.) 
Experience with wind tunnels has shown that the fiow will indeed self-adjust 
(in particular, the location of the acoustic horizon will self-adjust) so as to 
satisfy this fine tuning. We can now calculate the acceleration of the fiuid at 
the acoustic horizon by adopting L'Hopital's rule. 



-c%A'I,/Ah 
(c%-2a^- 



Now use the fact that 
Therefore 



-P 



{A v)' 
dp2 A V 

A" /A 



-PH 



d^p 
d^ 





an 


H 





2 + p{d^p/dp^)/c^ 



H 



(2.19) 



(2.20) 



(2.21) 
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That is 



aH 



2A \ 



A" 



(2.22) 



H 



Thus the physical acceleration of the fluid as it crosses an acoustic horizon 
is tightly constrained in terms of the speed of sound, the geometry of the 
horizon [A^ and A"^), plus some information coming from the equation of 
state. 



2.3.2 "Surface gravity" 



It is more useful to consider the "surface gravity" defined by the limit of the 
quantity 



' = --2 



ld(c2 



(2.23) 



It is this combination g, rather than the physical acceleration of the fluid a, 
that more closely tracks the general relativistic notion of "surface gravity", 
and it is the limit of this quantity as one approaches the acoustic horizon 
that enters into the Hawking radiation calculation fl^. Note that 



(2.24) 



This implies, in particular, that the fine-tuning (|2.18|) used to keep a finite 
at the acoustic horizon will also keep g finite there. In particular 



and so 



9h 



9h = clh 



1 + 



p d^p 
2?dp2 



H 



V2A 



H 




2c2 d/72 



(2.25) 



(2.26) 



H 



The first factor is of order / i?, with R the minimum radius of the nozzle, 
while the second and third factors are square roots of dimensionless numbers. 
This is in accord with our intuition based on dimensional analysis 0, [16| . If 
A" < 0, corresponding to a maximum of the cross section, then an and 
gn are imaginary which means no event horizon can form there. The two 
signs ± correspond to either speeding up and slowing down as you cross the 
horizon, both of these must occur at a minimum of the cross sectional area 
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A" > 0. (If the flow accelerates at the horizon this is a black hole horizon 
[future horizon]; if the flow decelerates there it is a white hole horizon [past 
horizon]. See Figure 1.) If the nozzle has a circular cross section, then the 
quantity A'^ is related to the longitudinal radius of curvature Rc at the throat 
of the nozzle, in fact 

A'i, = vr^. (2.27) 



2.3.3 Bose-Einstein condensate 

The technological advantages provided by the use of BECs as a working fluid 
for acoustic black holes have been discussed by Garay et al. (see also 



reference |T3| for a discussion of plausible extensions to that model). The 
present discussion can be interpreted as a somewhat different approach to 
the same physical problem, side-stepping the technical complications of the 
Bogoliubov equations in favour of a more fluid dynamical point of view. For 
a standard BEG 

m 

Then 



P 



P 



d(c^ 

dp 



while 



d^ 
d^ 



3 
2' 



So we have, rather simply 



(2.28) 



(2.29) 



(2.30) 



Similarly 



an 



9h 



H 



lA 



H 



3A^/4. 



(2.31) 



(2.32) 



This implies, at a black hole horizon [future horizon], a Hawking tempera- 
ture i, % ^ ^ 



k^TH 



h- 



ch 



l3A'i, 



27CCH l-K^/Au 



(2.33) 
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Ignoring the issue of gray-body factors (they are a refinement on the Hawking 
effect, not really an essential part of the physics), the phonon spectrum peaks 
at 

^peak = T^r4^ (2-34) 



that is 

Apeak = 47r2V^^/^. (2.35) 

This extremely simple result relates the typical wavelength of the Hawking 
emission to the physical size of the constriction and a factor depending on the 
flare-out at the narrowest point. Note that you cannot permit A'^j to become 
large, since then you would violate the quasi-one- dimensional approximation 
for the fluid flow that we have been using in this note. (There is of course 
nothing physically wrong with violating the quasi-one-dimensional approxi- 
mation, it just means the analysis becomes more complicated. In particular, 
if there is no external body force and the viscous forces are zero then by 
slightly adapting the analysis of [1^ the acoustic horizon [more precisely the 
ergo-surface] is a minimal surface of zero extrinsic curvature.) The preceding 
argument suggests strongly that the best we can realistically hope for is that 
the spectrum peaks at wavelength 

Apeak ~ \/^. (2.36) 

(Note that this is the analog, in the context of acoustic black holes, of the 
fact that the Hawking flux from general relativity black holes is expected 
to peak at wavelengths near the physical diameter of the black hole, its 
Schwarzschild radius — up to numerical factors depending on charge and 
angular momentum.) You can (in principle) try to adjust the equation of 
state to make the second factor in ( p.26| ) larger, but this is unlikely to be 
technologically feasible. 



2.3.4 Physical estimates 

It is the fact that the peak wavelength of the Hawking radiation is of order 
the physical dimensions of the system under consideration that makes the 
effect so difficult to detect. This suggests that it might be useful to look for 
indirect effects. In particular, in BECs it is common to have a sound speed of 
order 6 mm/s. If one then chooses a nozzle of diameter about 1 micron, and 



10 



a flare-out of A'^j ^ 1, tlien ~ 7 n K. Compare tliis to tlie condensation 
temperature required to form the BEC 

Tcondensate ^ 90 U K. (2.37) 

We see that in this situation the Hawking effect, although tiny, is at least 
comparable in magnitude to other relevant temperature scales. Moreover re- 
cent experiments indicate that it is likely that these figures can be improved. 
In particular, the scattering length for the condensate can be tuned by mak- 



ing use of the so called Feshbach resonance ||T8[. This effect can be used to 



increment the scattering length; factors of up to 100 have been experimen- 



tally obtained |jT9[- Therefore the acoustic propagation speed, which scales 
as the square root of the scattering length, could thereby be enhanced by a 
factor up to 10. This suggests that it might be experimentally possible to 
achieve ~ 6 cm/s, and so 

Th ^ 70 n K; (2.38) 

which places us much closer to the condensation temperature. The speed 
of sound can also be enhanced by increasing the density of the condensate 
(propagation speed scales as the square root of the density). In all of these 
situations there is a trade-off: For fixed nozzle geometry the Hawking tem- 
perature scales as the speed of sound, so larger sound speed gives a bigger 
effect but conversely makes it more difficult to set up the supersonic flow. 

The current analysis is purely "hydrodynamic" , and does not seek to deal 
with the "quantum potential" — the fact that the dispersion relation is at 
high momenta modifled in such a way as to recover "inflnite" propagation 
speed as in the Bogoliubov dispersion relation This issue has relevance to 
the so-called trans-Planckian problem (which in this BEC condensate context 
becomes a trans-Bohrian problem). Fortunately it is known, thanks to model 
calculations in fleld theories with explicit high-momentum cutoffs, that the 
low energy physics of the emitted radiation is largely insensitive to the nature 
and specific features of the cutoff. 

To summarize: this analysis complements that of Garay et al. [|T3|, in 
that it provides a rationale for simple physical estimates of the Hawking 
radiation temperature without having to solve the full Bogoliubov equations. 
Additionally, the current analysis provides simple numerical estimates of the 
size of the effect and identifies several specific physical mechanisms by which 
the Hawking temperature can be manipulated: via the speed of sound, the 
nozzle radius, the equation of state, and the degree of fiare-out at the throat. 
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2.4 Bogolubov dispersion relation 



However, there is a bit of a puzzle hiding in this analysis: We started with 
the nonlinear Schrodinger equation. That equation is parabolic, so we know 
that the characteristics move at infinite speed. How did we get a hyperbolic 
d'Alembertian equation with a finite propagation speed? The subtlety resides 
in neglecting the higher- derivative term Vq. To see this, keep Vq, and go to 
the eikonal approximation. One obtains the dispersion relation [l^, 



(uj -Vq - Pj 



cl e + 




(2.39) 



This is the curved-space generalization of the well-known Bogolubov disper- 
sion relation. Equivalently 



\ 



The group velocity is 



dcu 

Vg = — 
' dk 

while for the phase velocity 

UJ k 



vo + 



clk^+i — P 
\ 2m 



+ 2m? ^ 



elk- + 



(2.40) 



(2.41) 



{v^-k)k + \cl + 



Am? 



k. 



(2.42) 



Both group and phase velocities have the appropriate relativistic limit at 
low momentum, but then grow without bound at high momentum, leading 
to an infinite signal speed and the recovery of the parabolic nature of the 
differential equation at high momentum, {k ^ kc = m Csjh] equivalently in 
terms of the acoustic Compton wavelength \ <^\c = h/{rn Cg).) 
To investigate the situation a little more deeply, consider: 



u{k) 



ml + + 



' k- 
2m^ 



(2.43) 



(This is equivalent to the original Bogolubov dispersion relation. We have 
set the background flow to zero. In BEG condensates mo = but there 
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are other condensed matter systems where it need not be zero. Additionally 
c = h = 1 for simplicity.) Then at low momenta {k <^ mo) the dispersion 
relation is Newtonian 

uik)=mo + ^ + Oik^). (2.44) 

while at intermediate momenta (mo ^ <^ ^oo) it is (approximately) rel- 
ativistic. Perhaps surprisingly at large momenta {k ^ moo) the dispersion 
relation again takes on Newtonian form 

'^ik) = - + moo + 0(A;-2), (2.45) 

2moo 

and explicitly deviates from Lorentz symmetry. (Even more complicated de- 
viations from Lorentz symmetry are possible, see for example reference [pl|] .) 

The implication is this: If we consider a mode that far away from the 
horizon has a wave vector k that is well inside the "phonon" region of the 
dispersion relation, and then follow that mode back until it approaches the 
horizon, then near the horizon k diverges and the mode leaves the "phonon" 
region. It enters the "particle" region of the dispersion relation. This is the 
analog, in this particular condensed matter context, of the so-called trans- 
Planckian problem of general relativistic black hole physics. Fortunately it 
is now realised that the low-fc far-from-the-horizon physics of the Hawking 
effect is largely insensitive to the precise details of how the dispersion relation 
is modified by high-/c near-horizon physics. It is only part of the near-horizon 
physics, specifically the "surface gravity" that is really important in regards 
to the Hawking effect. 

As a closing comment we would like to add that dispersion relation of the 
form ( p. 43 ), and with quadratic or cubic deviations from Lorentz invariance. 



have also been encountered in several approaches to quantum gravity (see 
e.g. ||2^) and that there have been recent attempts to test these ideas via 



astrophysical observation. (See e.g. p3|, Ell and references therein. 



3 Slow light 

Slow light systems, photon pulses with anomalously low group velocities en- 
gendered by electromagnetically induced transparency (EIT) in an otherwise 
opaque medium, have also been mooted as being experimentally interesting 
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avenues towards building analogue black holes ^ One of the key 
issues here is that EIT intrinsically requires one to work in a narrow fre- 
quency range close to an atomic resonance; the resulting analogue black hole 
will trap pulses of light only over a very narrow frequency range, outside of 
which the medium is typically opaque. 

Although the technology for building and manipulating slow light systems 
is developing at an extremely rapid pace this intrinsic limitation to 
working in a narrow frequency range somewhat obscures the meaning of 
Hawking radiation and makes it less clear just what signal should be looked 
for. For a discussion of the possibilities see p8[1 . 



4 Quasiparticles 

The use of superfluid quasiparticles, in particular the quasiparticles and do- 



main walls of liquid HesA, has been investigated by Jacobson and Volovik |29 
A particularly nice feature is due to the two-fluid nature of the system, in 
that in this system it seems possible to arrange a wide separation between 
the Landau critical velocity and the velocity relevant to defining the horizon 
at which the Hawking phenomenon is expected to occur. 



5 Normal modes 

The sheer number of different physical systems in which analogue models for 
general relativity may be found is indicative of a deep underlying principle. 
Indeed, finding an approximate Lorentzian geometry is really just a matter of 
picking an arbitrary physical system, isolating a particular degree of freedom 
that is approximately decoupled from the rest of the physics, and doing a 
low-momentum field-theory normal-modes analysis [^, 

Roughly speaking: in any hyperbolic system of differential equations (no 
matter how derived) there are by definition wave-like solutions P2|, 53 . The 



set of admissible wavevectors associated with these wave-like solutions can 
be used to define (modulo some nasty complications we defer to the tech- 
nical literature |Q) a cone-like structure in momentum space, and hence a 
conformal class of Lorentzian-signature metrics. For this reason the emer- 
gence of Lorentzian-signature effective metrics is an almost generic aspect of 
low-momentum physics. 
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6 Emergent gravity 



So far, the entire discussion has been about models of gravity, models that 
reproduce the kinematics . If one wants to make a bolder proposal, that 
analogue models might be useful for generating models for gravity, mod- 
els that reproduce the Einstein-Hilbert dynamics (or some approximation 
thereto), then the situation is considerably more subtle and tentative. Kine- 
matics is relatively easy, and is in some sense generic. Einstein-Hilbert dy- 
namics is trickier — to get an "emergent gravity" arising from these ana- 
logue models will require some variant of Sakharov's notion of "induced 
gravity" A useful observation in this regard is that any curved-space 

relativistic quantum field theory will automatically generate an Einstein- 
Hilbert counterterm through one- loop effects In heat kernel language, 
the first Seeley-DeWitt coefficient generically contains a term proportional 
to the Einstein-Hilbert action, and after renormalization this generically pro- 
vides an Einstein-Hilbert term in the effective action |^ . Unfortunately the 
same logic provides an uncontrolled cosmological constant from the zeroth 
Seeley-DeWitt coefficient, plus quadratic curvature-squared terms from the 
second Seeley-DeWitt coefficient, so the argument is not fully acceptable. 
Furthermore there are technical issues involved in specifying the volume of 
the function space on which this effective action is defined. To get Ein- 
stein gravity one needs both an Einstein-Hilbert action and the freedom to 
perform arbitrary metric variations. Though the situation is still far from 
clear, interest in these possibilities is both long-standing (see for instance the 
sub- manifold models in references pTl and ongoing Bl, |3^, HO. 



7 Discussion 

In this mini-survey we have seen how an effective metric emerges as a low- 
energy low-momentum approximation in certain physical systems. Indeed 
we have been able to argue that the emergence of such effective metrics is an 
almost generic consequence of performing a "normal modes" analysis on an 
arbitrary field theory. Once one has an effective metric in hand (no matter 
how derived), all kinematic aspects of general relativity can in principle be 
carried over to these analogue systems — in particular all curved-space field 
theory (both classical and quantum) finds a natural home in these systems. 
The most stunning feature of these analogue models is the ability to 
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generate analogue horizons (analogue black holes) and more specifically, the 
possibility of detecting an analogue form of Hawking radiation. In the body 
of this article we have specifically considered the use of acoustics in Bose- 
Einstein condensates as a particularly promising analogue system. This par- 
ticular model stands out for purely technological reasons — the condensation 
temperature, required to form the condensate in the first place is of order 
90 nK, which is considerably less than 1 order of magnitude away from the 
estimates of the relevant Hawking temperature. It is this congruence between 
two important physical scales that makes this particular system so interest- 
ing. Many condensed matter systems are capable of mimicking curved space 
quantum field theory; this particular condensed matter system does so in a 
particularly interesting manner that seems amenable to experimental probes 
in the not too distant future. 
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